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TECHNICAL  MEMORANDUM 


NIGHTTIME  OBSERVATIONS  OF  THUNDERSTORM  ELECTRICAL  ACTIVITY 

FROM  A  HIGH  ALTITUDE  AIRPLANE 

INTRODUCTION 


Airplane  and  balloon  observations  have  confirmed  C.  T.  R.  Wilson's  prediction 
that  strong  electric  fields  exist  over  the  tops  of  thunderstorms  [2].  Photographs 
from  airplanes  made  in  the  daytime  in  conjunction  with  electric  field  measurements 
show  that  the  most  intense  electrical  stresses  are  associated  with  the  penetrative 
convective  cells  that  rise  above  the  cirriform  anvil  [3,4]. 

At  the  tops  of  thunderstorms  there  is  undoubtedly  strong  electrical  activity, 
but  there  are  only  a  few  observations  concerning  lightning  or  other  kinds  of  dielectric 
breakdown  that  may  be  taking  place.  Descriptions  of  discharges  in  the  clear  air 
above  the  clouds  are  to  be  found  in  texts  on  atmospheric  electricity.  Schonland  [5] 
tells  of  a  brilliant  ball  of  light  observed  at  the  top  of  a  bank  of  clouds  which  pulsated 
in  size  and  intensity  for  over  15  min.  He  also  described  what  appears  to  be  a  form 
of  glow  discharge  which  extended  from  the  top  of  a  cloud  to  the  upper  air  and  was 
seen  to  occur  regularly  at  the  same  time  as  flashes  to  ground. 

Chalmers  [6]  cites  visual  evidence  of  possible  discharges  toward  the  electro¬ 
sphere  reported  by  Boys,  Malan,  Wood,  Reynolds,  and  Hoffman,  as  well  as  observa¬ 
tions  of  radar  reflections  by  Rumi  and  Atlas,  which  may  have  originated  from  such 
discharges.  Uman  [7]  cites  three  additional  observations  of  discharges  that  extended 
upward  from  the  tops  of  thunderclouds. 

In  private  communications  we  have  learned  of  similar  independent  observations 
of  glow-like  discharges  above  severe  storms  from  C.  B.  Moore  of  New  Mexico  Tech 
and  from  Floyd  Montgomery,  formerly  an  observer  for  the  U.S.  Weather  Bureau. 
Recently,  Ronald  Williams  [8]  has  described  a  lightning-like  discharge  from  the  top 
of  a  cloud  into  the  stratosphere  which  he  saw  at  night  while  flying  a  U-2  airplane  at 
approximately  20  km  altitude.  Similar  observations  from  airplanes  by  Roberts  [9]  and 
Gales  [10]  have  been  published. 

Possibly  related  to  such  observations  is  a  phenomenon  sometimes  observed  in 
which  the  clear  atmosphere  above  a  vigorous  storm  can  be  seen  to  give  off  a  flickering 
light.  This  may  possibly  be  light  from  the  lightning  discharges  that  is  being 
scattered  by  aerosol  particles  in  the  stratosphere.  Alternatively,  it  is  conceivable 
that  such  light  may  originate  as  the  result  of  glow-like  electrical  discharges  in  the 
rarified  upper  atmosphere.  That  such  discharges  could  exist  in  the  clear  air  above 
a  thundercloud  seems  plausible,  for  the  largest  and  most  electrically  active  thunder¬ 
clouds  are  known  to  rise  to  altitudes  in  excess  of  20  km ,  where  the  atmospheric 
pressure  and  hence  the  dielectric  breakdown  strength  is  less  than  l/10th  that  at  sea 
level.  Furthermore,  as  C.  T.  R.  Wilson  [11]  has  pointed  out,  the  dielectric  strength 
of  the  air  will  decrease  more  rapidly  with  altitude  than  will  the  electric  field  perturba¬ 
tions  that  result  from  lightning.  It  therefore  appears  conceivable  that  discharges 
might  take  place  in  the  stratosphere  some  distance  above  the  cloud. 


Although  the  number  of  observations  of  discharges  in  the  clear  air  above  the 
cloud  is  small,  they  probably  occur  more  often  than  might  at  first  be  supposed,  for 
observers  rarely  have  a  clear  view  of  a  thunderstorm  cloud  top.  Schonland  [5]  may 
be  correct  in  his  suggestion  that  though  this  type  of  air  discharge  is  rarely  reported, 
it  may  not  be  uncommon. 

Photographs  of  the  lightning  associated  with  the  upper  part  of  a  thundercloud 
would  be  of  help  in  understanding  this  phenomenon,  but  they  are  rare.  For  example, 
no  such  picture  of  lightning  is  to  be  found  in  Salanave's  beautiful  new  atlas  of 
lightning  photographs  [12].  The  dearth  of  such  pictures  is  not  surprising.  The 
upper  parts  of  large  storms  are  usually  shielded  from  the  view  of  ground  observers 
by  low-level  clouds  or  by  widespreading  anvils.  From  aircraft,  photography  looking 
down  at  cloud  tops  is  usually  not  feasible,  for  the  most  intense  electrical  storms 
attain  heights  well  above  the  ceilings  of  all  but  a  few  specialized  high-altitude  air¬ 
planes.  Even  when  the  top  of  the  cloud  is  in  clear  view,  either  from  the  ground  or 
from  an  airplane,  photographing  lightning  during  the  day  poses  problems,  for  the 
lightning  channel  is  hardly  visible  when  viewed  against  the  sky  or  the  bright  back¬ 
ground  of  a  sunlit  cloud. 

There  can  be  little  question  that  a  better  knowledge  of  the  electrical  activity  in 
the  upper  part  of  the  cloud  and  its  relationship  to  the  convective  structure  of  the 
storm  is  important  to  an  understanding  of  thunderstorms.  These  discharges  may  not 
only  be  important  in  the  electrical  budget  of  the  cloud,  but  may  also  play  a  significant 
role  in  introducing  oxides  of  nitrogen  and  other  chemical  compounds  directly  into  the 
stratosphere. 

Astronauts  in  the  past  have  described  spectacular  lightning  occurring  in  the 
tops  of  thunderstorms  [  13] .  These  observations  show  that  an  orbiting  spacecraft  is 
a  valuable  platform  from  which  to  observe  the  electrical  activity  in  clouds.  To  pro¬ 
vide  more  detailed  data  from  such  observations,  we  have  made  arrangements  to  have 
astronauts  on  space  shuttle  flights  take  pictures  of  thunderstorms  and  lightning  with 
a  motion  picture  camera  and  to  utilize  a  photocell  optical  system  to  secure  information 
on  the  frequency  and  character  of  lightning  discharges  [14,15].  More  recently,  the 
possibility  of  making  meaningful  lightning  measurements  from  a  geosynchronous 
satellite  has  been  discussed  [  16] . 

In  order  to  provide  background  data  and  ground  truth  information  that  would 
be  helpful  in  the  interpretation  of  satellite  data  and  in  the  design  of  new  satellite 
instrumentation,  we  have  made  several  flights  over  active  thunderstorm  cells  with  a 
NASA  U-2  airplane  instrumented  for  electrical  and  optical  observations.  The  results 
of  the  first  of  these  flights,  which  was  made  over  a  thunderstorm  in  New  Mexico 
during  the  day,  have  been  reported  [17]. 

This  report  describes  two  subsequent  sets  of  observations  that  were  made 
looking  down  at  night  from  the  U-2  airplane,  which  was  flying  at  an  altitude  of  about 

20  km  and  a  true  air  speed  of  200  m  sec  The  first  of  these  on  17  May  1980  was 
made  over  a  vigorous  thunderstorm  taking  place  160  km  northwest  of  Little  Rock, 
Arkansas,  and  the  second  on  11  May  1981  over  a  thunderstorm  near  Atlanta,  Georgia. 
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INSTRUMENTATION 


Cameras 

Three  Vinten  cameras,  normally  employed  for  aerial  photography  of  the  ground 
from  the  U-2  aircraft,  were  used  to  photograph  thunderstorms  and  lightning.  Each 
of  these,  which  looks  directly  down  from  the  airplane,  has  an  f/2.8  lens  of  44.5  mm 
focal  length  and  uses  70  mm  film,  giving  a  format  of  56  x  57  mm.  When  the  airplane 
is  at  its  operating  altitude  of  20  km,  the  ground  coverage  is  about  26  x  26  km.  For 
clouds  at  intermediate  altitudes,  such  as  are  described  in  this  report,  the  coverage 
is  approxiamtely  10  x  10  km. 

For  photography  of  the  Arkansas  storm,  all  cameras  were  loaded  with  black  and 
white  film ,  the  first  with  Plus-X  Aerographic  2402  to  bring  out  cloud  details ,  the 
second  with  Infrared  Aerographic  2424  to  determine  if  this  film  offered  any  particular 
advantages  for  cloud  photography  at  night,  and  the  third  with  Panatomic-X  Aerial 
3400  film  to  give  the  maximum  resolution  of  the  lightning  discharge.  This  camera  was 
also  equipped  with  a  300  lines /mm  diffraction  grating  placed  over  the  lens  to  provide 
slitless  spectra  when  lightning  channels  were  visible  in  the  field  of  view. 

For  the  daytime  photography  of  thunderstorms,  the  cameras  were  operated  in 
their  normal  mode,  exposing  snapshots  in  unison,  and  advancing  the  film  between 
exposures.  For  night  photography  of  thunderstorms  and  lightning  the  shutters  were 
removed .  In  the  flight  over  the  Arkansas  storm ,  the  cameras  were  actuated  by  the 
pilot,  a  procedure  which  proved  satisfactory  for  lightning  photography. 

The  photographic  procedures  were  modified  slightly  for  the  flight  over  the 
second  storm  made  near  Atlanta,  Georgia.  An  intervalometer  was  used  to  expose 
each  frame  for  11  sec.  In  one  of  the  cameras,  Eastman  SO 397  color  film  was  used 
instead  of  high-speed  black  and  white  film  to  facilitate  the  interpretation  of  spectra 
produced  by  the  diffraction  gratings,  which,  on  this  flight,  were  placed  over  the 
lenses  of  the  color  camera  and  of  the  infrared  camera. 


Photocell  Transient  Detector 

A  silicon  solar  cell  looking  downward  through  a  glass  window  was  used  to  pro¬ 
vide  electrical  signals  for  recording  the  details  of  the  optical  transients  produced  by 
lightning.  A  red  to  near  infrared  optical  filter  (Baird  Atomic  D-645)  having  85  per¬ 
cent  of  maximum  transmission  at  X  =  656  nm  was  placed  in  front  of  the  photocell. 

An  absolute  calibration  of  the  filter  plus  photocell  unit  gave  a  sensitivity  of  0.31  A/W 
into  a  5  ohm  load.  This  equipment  is  identical  to  that  which  was  used  in  daylight 
observations  of  lightning  from  a  U-2  [17].  High  gain  and  low  gain  channels  were 
used  to  record  the  large  dynamic  range  of  the  signals. 


Electric  Field-Change  Meter 

Because  measurements  of  the  electric  field  changes  are  valuable  in  the  diagnosis 
of  the  various  processes  taking  place  in  the  lightning  flash ,  an  electric  field-change 
meter  was  also  installed  behind  the  down-facing  window  of  the  U-2  aircraft.  The 

o 

meter  consisted  of  a  flat-plate  sensor  of  area  300  cm  connected  to  a  charge  amplifier 
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with  a  3  sec  time  constant.  Three  output  channels  set  to  different  gains  were  used 
to  cover  the  dynamic  range  of  the  E-field  signals.  The  rise-time  of  the  meter  from 
10  to  90  percent  amplitude  was  less  than  100  nanosec. 


Linear  Diode  Array  Spectrometer 

In  order  to  supplement  the  spectroscopic  data  being  recorded  on  film,  a  line- 
scan  camera  employing  a  300  lines /mm  grating  and  a  1728  linear  diode  array  was 
installed  in  the  lower  bulkhead.  This  equipment  (Fairchild  Line-Scan  Camera  Model 
CCD  1300)  generated  spectral  data  as  a  video  signal  that  was  recorded  along  with  the 
photocell  and  the  E-field  data  on  a  Bell  and  Howell  M-14E  tape  recorder.  Data  were 
taken  at  a  tape  speed  of  30  ips.  The  slow  tape  speed  compromised  the  bandwidth 
available  with  the  wideband  II  electronics,  but  was  used  in  order  to  achieve  longer 
total  recording  time. 

The  optical  and  the  E-field  data  recorded  at  30  ips  do  not  exhibit  the  lightning 
fine-structure  present,  for  example,  in  the  E-field  data  of  Guo  and  Krider  [18]  who, 
using  ground  based  equipment ,  studied  the  relationship  of  the  light  emitted  to  the 
accompanying  electric  field  signatures  of  return-strokes.  In  general,  it  was  found 
that  the  rise  time  of  the  optical  signals  was  slow ,  much  slower  than  the  measured 
10  to  90  percent  value  of  1.6  nsec  of  the  optical  sensor,  and  comparable  with  the 
rise  times  associated  with  the  return  strokes  as  given  in  Reference  18.  There  are 
several  reasons  for  the  apparent  dilation  of  the  optical  signals  when  measured  with 
a  sensor  looking  down  at  the  thundercloud:  the  leader  and  the  return-stroke 
channels  are  usually  not  seen  directly,  but,  as  they  extend  their  length,  multiple 
scattering  coupled  with  the  large  dimension  of  the  cloud  and  the  finite  velocity  of 
light  results  in  a  lengthening  of  the  light  emitting  channel  and  in  multipaths  that 
contribute  to  a  "stretched  signal."  Thus,  lightning  optical  studies  that  use  a 
horizontal  slit  to  view  a  small  portion  of  the  luminous  channel,  such  as  discussed  in 
Beasely,  et  al.  [19]  probably  measure  the  true  optical  risetime  in  that  portion  of  the 
channel.  The  wide  angle  optical  detector  (120  deg)  used  on  the  U-2  viewing 
primarily  scattered  light  emerging  from  the  cloud  top  was  probably  more  of  a  deter¬ 
mining  factor  in  the  measured  slow  rise  times  than  was  the  reduced  recording  speed. 
On  the  other  hand,  since  one  major  objective  of  this  study  was  simply  to  resolve 
leaders  from  return- strokes ,  the  30  ips  tape  speed  was  more  than  adequate  for  the 
purpose.  The  lack  of  greater  detail  in  the  E-field  signal  is  due  in  part  to  the 
necessity  for  smoothing  the  data  to  reduce  the  large  amount  of  spurious  noise  picked 
up  by  the  meter  from  various  sources  in  the  U-2  aircraft  bay. 


OBSERVATIONS 


The  Arkansas  Storm 

At  0000  GMT  17  May  1980  (about  3  hr  prior  to  flight  time),  the  surface  weather 
map  featured  a  1006  mb  low  centered  in  south  central  Kansas.  A  trailing  cold  front 
extended  southeastward  across  extreme  eastern  Oklahoma  and  then  southward  into 
Texas.  A  warm  front  extended  eastward  from  the  low  across  southern  Missouri. 

The  flight  target  area  of  the  U-2  aircraft  was  in  the  apex  of  the  warm  section  of  the 
storm.  Farther  south  an  intense  mesoscale  convective  complex  was  located  along  the 
southern  Louisiana  coast  and  was  responsible  for  4  to  8  cm  rains  in  the  last  6  hr  at 
several  locations. 
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Surface  winds  in  the  target  area  were  southeasterly  between  2  to  4  m  s  at 
0000  GMT.  Surface  dewpoints  were  near  20°C  with  air  temperatures  around  25°C. 

A  local  moisture  maximum  was  observed  at  850  mb  at  Little  Rock  (72340)  with  a  mixing 

ratio  of  12  g  kg  ^  and  a  southerly  wind  of  10  m  s  In  the  upper  troposphere  a 
broad  southwesterly  flow  was  evident  ahead  of  a  vigorous  trough  located  in  western 

Texas.  200  mb  winds  were  southwesterly  at  45  m  s  *  over  the  target  area. 

The  Arkansas  convection  formed  just  ahead  of  the  cold  front  and  on  the  back 
edge  of  the  intense  convective  complex  to  the  southeast.  The  Showalter  index  at 
0000  GMT  17  May  1982  was  0,  -5  and  -5  at  Longview,  Texas  (72247),  Little  Rock, 
Arkansas  (72340),  and  Monett,  Missouri  (72349),  respectively.  The  most  unstable  air 
was  found  over  western  Arkansas  and  southwestern  Missouri  where  the  low  level 
moisture  maximum  was  capped  by  dry  air  in  the  700  to  500  mb  layer. 


Photographic  Data 

On  this  flight,  the  camera  was  operated  by  the  pilot.  When  he  flew  over  the 
storm  and  if  lightning  was  occurring,  he  would  actuate  the  camera.  Because  the 
LED  timer  was  not  operative,  the  time  of  each  exposure  is  not  known,  and  correla¬ 
tions  cannot  be  made  between  the  photographs  and  the  data  taken  with  the  other 
instruments. 

When  the  lightning  within  the  storm  was  sufficiently  intense,  the  entire  cloud 
surface  in  the  field  of  view  of  the  camera  was  illuminated  and  revealed  the  cauliflower¬ 
like  structure  typical  of  a  strongly  convective  turret.  Figure  1  shows  such  a  cloud. 
The  distance  scale  given  is  only  approximate,  as  it  is  calculated  from  the  height  of 
the  cloud ,  which  was  estimated  by  the  pilot  to  be  at  about  14  km . 

In  most  photographs,  the  lightning  channels  lie  entirely  beneath  the  visible 
cloud  surface,  and  cannot  be  seen.  In  some  photographs,  perhaps  10  percent  of  the 
total,  a  portion  of  the  sharp  bright  channel  of  an  electric  discharge  is  clearly  visible. 
Figure  2  shows  the  most  common  situation  in  which  a  section  of  the  channel  can  be 
seen  coming  out  of  one  part  of  the  cloud,  crossing  an  open  crease  or  valley  in  the 
cloud,  and  then  re-entering  the  cloud  on  the  opposite  side.  Two  strokes  in  the  same 
channel  are  visible  because  the  discharges  occurred  while  the  film  was  being  trans¬ 
ported  past  the  unshuttered  lens.  Occasionally,  as  is  shown  in  Figure  3,  a  lightning 
channel  can  be  seen  to  terminate  in  the  clear  air.  Figure  4  shows  a  bright  spot, 
which  is  perhaps  produced  when  the  light  escapes  through  a  hole  in  the  cloud. 


Optical  and  Electric  Field-Change  Data 

Figure  5(0)  shows  the  optical  and  electric  field  signal  for  a  lightning  flash  with 
5  strokes  to  ground,  which  occurred  on  Day  138  at  02:25:22  in  the  morning.  Because 
we  have  no  ground-truth  data,  the  identification  of  the  flash  as  a  flash  to  ground  is 
an  interpretation  based  on  the  similarity  of  these  field  changes  to  others  recorded 
and  verified  at  ground  level.  The  E- field-change  record  shown  in  Figure  5(0) 
exhibits  5  downward- going  abrupt  deflections  each  indicative  of  a  return-stroke 
equivalently  lowering  negative  charge  to  Earth  (or  raising  positive  charge  to  the 
cloud).  Each  stroke  is  labeled  in  the  order  of  occurrence,  and  is  shown  on  an 
expanded  time  base  in  Figures  5(1)  to  5(5).  The  abrupt  E-field  changes  are  each 
accompanied  by  an  optical  signal,  which  is  also  shown  expanded  in  the  figures. 
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Figure  1. 
15  km 


Convective  cloud  turret  of  Arkansas  Storm  (approximately 
msl  altitude)  illuminated  by  lightning  within  the  cloud. 
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Figure  2.  Lightning  channel  in  clear  air  bridging  crease  in  cauliflower  turret. 
Double  image  is  produced  by  two  strokes  in  the  same  channel  that 
took  place  during  transport  of  film  past  the  open  shutter. 
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Figure  3.  Lightning  channels  apparently  terminating  in  clear  air. 
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Figure  4.  Bright  spot  of  light  in  dark  cloud. 
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Figure  5.  Optical  sensor  and  electric  field-change  data  from  a  flash  to  ground 
near  Little  Rock,  Arkansas.  The  U-2  aircraft  was  flying  at  an  elevation  of 
20  km.  Figure  5(0)  shows  the  sensor  data  for  the  entire  flash. 

Figure  5(1)  through  5(5)  are  expanded  views  of  the 
signals  from  the  individual  strokes. 
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Focusing  attention  on  the  first  stroke,  in  Figure  5(0),  it  is  noted  that  the 
E-field  signal  shows  a  small  downward-going  trend  for  about  15  msec  prior  to  the 
return  stroke.  This  deflection  is  probably  associated  with  a  stepped  leader:  not 
only  does  it  appear  before  the  first  stroke  of  the  flash,  but  its  duration  of  15  msec 
is  very  much  what  is  expected  for  stepped  leaders.  In  the  expanded  time  base  of 
Figure  5(1),  about  1.5  msec  of  the  leader  field  change  is  shown.  About  600  ysec  of 
activity  is  discernible  in  the  E-field  record  immediately  prior  to  the  stroke,  and  this 
appears  primarily  as  a  smooth  downward  trend  without  pulsative  noise  because,  as 
explained  earlier,  the  noise  environment  required  that  the  E-field  data  be  smoothed. 
The  optical  signal  accompanying  the  stepped  leader  was  apparently  too  weak  to  be 
recorded.  The  optical  signal  for  the  return-stroke  portion  is  strong,  and  exhibits  a 
complex  structure  for  a  duration  somewhat  less  than  1  msec.  This  type  of  structure 
was  attributed  by  Guo  and  Krider  [18]  to  the  occurrence  of  branches.  From  the 
above-cloud  view,  it  is  difficult  to  attribute,  with  certainty,  the  obviously  multiple 
and  superimposed  signals  to  a  specific  source  such  as  below-cloud  branching  in  the 
return-stroke  channel.  As  seen  from  above  the  cloud,  below-cloud  branching  would 
be  very  strongly  attenuated.  In-cloud  branching  is  perhaps  a  more  likely  source. 

In  Figure  5(2),  the  second  return- stroke  is  preceded  by  a  leader  of  about 
1.6  msec  duration  as  seen  in  the  E-field  trend.  The  optical  return-stroke  signal  is 
also  preceded  by  a  light  pulse  that  occurs  simultaneously  with  the  E-field  leader 
signal.  The  1.7  msec  separation  between  the  leader  and  return-stroke  optical  signals 
is  somewhat  long  but  still  consistent  with  the  well  established  dart-leader,  return- 
stroke  sequence.  The  dart  leader  optical  signal  peak  and  the  following  return  to 
zero  —  all  prior  to  the  return-stroke  —  is  interesting  and  will  be  discussed  presently. 

The  third  stroke,  Figure  5(3),  is  an  excellent  example  of  the  dart-leader, 
return-stroke  sequence  in  both  the  optical  and  E-field  signals.  The  large  optical 
return-stroke  signal  is  seen  to  go  off-scale.  In  this  sequence,  the  dart  leader  peak 
occurred  about  500  ysec  prior  to  the  return-stroke  peak. 

The  fourth  stroke  exhibits  a  complex  optical  return-stroke  signal  preceded  by 
a  discernible  luminosity.  There  seems  to  be  no  accompanying  E-field  deflection  prior 
to  the  return- stroke;  the  long  time  interval  between  the  third  and  fourth  stroke 
(120  msec)  and  the  sporadic  but  low  level  optical  activity  in  that  interval  is  suggestive 
of  a  second  stepped-leader .  In  this  case,  again,  the  complex  structure  of  the  return- 
stroke  optical  signal  could  be  representative  of  strong  branching  associated  with  a 
first  stroke  in  a  new  channel. 

In  Figure  5(5)  a  clear  example  of  an  optical  dart-leader,  return-stroke  sequence 
is  seen. 

In  the  second ,  third ,  and  fifth  strokes  light  signals  are  identified  as  dart 
leaders  because  (1)  their  luminous  intensity  is  less  than  the  intensity  of  the  return- 
stroke  that  follows,  (2)  they  occur  about  1.5  msec  or  less  prior  to  the  return-stroke, 
and  (3)  they  correlate  well  with  the  leader  E-field  change.  Nevertheless,  the  optical 
signals  associated  with  the  dart  leader  all  reach  a  peak  and  then  decrease  before  the 
return-stroke  luminosity  appears  as  a  well  developed,  fast-rise  signal.  This  is  not, 
at  first  sight,  what  one  might  expect  a  dart  leader  optical  signal  to  look  like.  There 
is  little  reason  to  expect  the  luminosity  to  decrease  so  abruptly  just  prior  to  the 
return- stroke;  in  fact,  one  might  expect  the  light  signal  to  increase  as  the  leader 
approaches  close  to  Earth.  We  believe  that  the  large  decrease  in  leader  luminosity 
is  not  due  to  any  change  in  the  emission  associated  with  the  leader ,  but  is  caused 
by  a  change  in  the  light  that  emerges  from  the  cloud  top  as  the  leader  traveling 
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downward  through  the  cloud  emerges  from  the  cloud  base.  As  shown  by  Thomason 
and  Krider  [1],  the  fraction  of  light  that  escapes  through  the  cloud  surface  is  high 
when  the  source  is  within  the  cloud,  the  primary  loss  being  due  to  absorption,  which 
is  negligible  in  the  visible  and  low  in  the  near  infrared.  On  the  other  hand,  the 
same  authors  show  that  for  a  light  source  below  a  cloud,  even  a  thin  cloud,  almost 
all  the  light  will  be  reflected  and  very  little  will  emerge  from  the  other  side.  With 
this  interpretation,  the  decrease  from  peak  luminosity  reached  by  the  leader  as 
indicative  of  the  emergence  of  the  leader  into  the  clear  air  below  the  cloud  is 
explained.  As  observed  from  scattered  light  by  the  U-2,  the  results  appear  to  be 
wholly  consistent  with  the  conclusions  reached  by  Thomason  and  Krider. 

The  E-field  change  and  the  optical  signals  from  an  intracloud  flash  that 
occurred  in  the  same  storm  at  02:31:05  are  shown  in  Figure  6.  This  flash  was  one 
of  short  duration  lasting  only  about  250  msec.  It  is  very  different  from  the  flash  to 
ground  shown  in  Figure  5(0)  in  the  pulsative  structure  of  the  optical  signal.  The 
upward  going  quasi-steady  deflection  in  the  E-field  also  differentiates  it  from  the 
ground  flash.  (The  noise  on  the  E-field  record  is  particularly  obvious  as  a  cyclic 
recurrence,  masking  any  small  abrupt  field  changes  that  may  have  been  present). 

An  additional  4  intracloud  flashes  and  14  cloud- to-ground  flashes  were  analyzed. 
Although  there  are  some  variations,  the  two  flashes  discussed  above  are  representa¬ 
tive  of  the  data  acquired  during  the  5/18/80,  Day  138,  Arkansas  flight. 
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Figure  6.  Optical  and  E  field-change  signal  from  an  intracloud  lightning 
flash  recorded  by  a  U-2  aircraft  flying  at  an  elevation  of  20  km 
above  a  storm  near  Little  Rock,  Arkansas. 
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Spectral  Data 


Results  of  the  photographic  spectral  data  are  summarized  in  Figures  7  and  8 
and  Table  1.  By  chance,  the  film  was  moving  at  the  time  of  the  flash,  so  that  the 
flash  was  resolved  into  two  channels:  one  with  five  strokes  and  the  other  with  four 
strokes  (Fig.  7).  Figure  8  shows  the  slitless  spectral  data  for  one  stroke  viewed 
with  a  video  densitometer  [20].  The  cursor  is  aligned  perpendicular  to  the  spectral 
lines,  and  the  resulting  density  profile  is  shown  below  the  cursor.  Note  carefully 
that  the  density  profile  is  not  calibrated  so  that  the  relative  intensity  of  the  lines 
cannot  be  expressed  quantitatively.  Nevertheless,  the  relative  brightness  can  be 
estimated. 

Table  1  summarizes  the  line  wavelengths,  their  qualitative  relative  intensity, 
and  their  identification.  In  a  few  cases  the  weak  lines  cannot  be  identified.  The 
most  intense  emissions  occur  at  463.0  nm  and  500.5  nm  and  are  from  the  confluence 
of  lines  from  singly  ionized  nitrogen  [21].  The  only  identified  neutral  emission  is 
from  H-beta  at  486.1  nm.  Reference  to  the  time-resolved  spectrum  of  Orville  (Ref. 
22,  Part  I,  Figure  4a)  suggests  that  the  lightning  channel  cooled  very  rapidly. 

This  is  suggested  on  the  basis  of  the  presence  of  singly  ionized  lines  and  the  lack 
of  significant  neutral  emission  lines.  The  latter  are  emitted  at  lower  temperatures 
and  are  characteristic  of  continuing  currents  [23].  Unfortunately,  it  is  not  known 
if  the  observed  flash  is  cloud-to- ground  or  intracloud. 

As  might  be  expected  from  the  spectral  data,  photographs  taken  with  the  color 
film  show  the  lightning  channels  to  be  blue  and  the  clouds  illuminated  by  the  lightn¬ 
ing  to  be  white  or  bluish  white. 

TABLE  1.  SPECTRAL  LINE  IDENTIFICATION 


Wavelength 

(nm) 


Probable 

Strength  Identification  and  Wavelength 
(estimated)  (nm) 


452.7 

m 

463.0 

s 

478.7 

m 

485.0 

m 

500.0 

s 

505.5 

m 

518.9 

s 

531.5 

w 

539.4 

w 

N 

II 

(59) 

453.0 

N 

II 

(5) 

463.0 

N 

II 

(20) 

480.3 

H- 

-beta 

486.1 

N 

II 

(19) 

500.5 

N 

II 

(4) 

504.5 

N 

II 

(66) 

517.9 

Line  strengths  are  listed  as  strong  (s),  medium  (m),  or 
weak  (w).  Probable  identifications  are  assigned  a  wave¬ 
length  corresponding  to  the  strongest  line  in  the  multiplet 
(Moore  [27]). 
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Figure  7.  Flash  occurring  while  film  was  moving  so  that  five  strokes  in  one  channel  and  four  strokes 
in  a  second  channel  of  the  flash  are  resolved.  Spectrum  from  the 
brightest  stroke  appears  in  a  corner. 


Figure  8.  Original  slitless  spectrum  negative  is  shown  on  a  video  densitometer. 
Density  profile  is  not  calibrated,  so  the  relative  changes 
show  only  qualitative  differences. 
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The  Georgia  Storm 


Three  hours  before  the  pictures  and  data  were  taken  over  the  thunderstorm 
near  Atlanta,  Georgia,  the  map  showed  a  low  in  southern  Indiana  (1,000  mb)  with  a 
cold  front  extending  south  across  southern  Indiana,  central  Kentucky,  central 
Tennessee,  eastern  Alabama  to  the  Florida  panhandle.  There  was  a  stationary  front 
running  east-northeast  across  southern  Indiana,  northeastern  Ohio,  northwestern 
Pennsylvania,  central  New  York,  northern  Vermont,  New  Hampshire  and  Maine.  A 
squall  line  existed  ahead  of  the  cold  front  in  extreme  eastern  Tennessee,  extending 
to  southwestern  North  Carolina  and  western  Georgia.  The  squall  line  was  oriented 
north  to  south  parallel  to  the  front  and  200  km  in  advance  of  it.  500  mb  winds  of 

30  m  sec  ^  out  of  the  south  southwest  prevailed  over  the  cold  front  and  squall  line. 
Representative  inflow  air  into  the  squall  line  had  a  surface  based  lifted  index  of  -6, 
indicative  of  the  ambient  unstable  conditions. 


Photographic  Data 

The  photographs  obtained  looking  down  on  this  storm  over  Georgia  showed 
features  similar  to  those  observed  over  Arkansas.  The  camera  was  again  operated  in 
the  open-shutter  model,  but  in  this  case  the  film  was  automatically  transported  to 
expose  each  frame  for  11  sec.  The  time  of  each  exposure  appears  on  the  frame,  and 
therefore  it  is  possible  to  correlate  the  photographs  with  other  data  that  were  taken. 
One  such  photograph  (Fig.  13),  will  be  discussed  in  connection  with  the  analysis  of 
the  optical  and  E-field  events  shown  in  Figure  12.  From  the  known  altitude  of  the 
airplane  and  air  speed,  it  is  possible  to  calculate  the  altitude  of  the  cloud  features 
by  their  apparent  motion  between  flashes  of  lightning.  In  this  storm,  such  calcula¬ 
tions  show  that  the  cloud  top  was  about  9  km  beneath  the  airplane,  or  at  an  altitude 
of  about  11  km  msl.  The  scale  indicated  on  these  photographs  is  probably  somewhat 
more  accurate  than  that  on  the  photographs  of  the  Arkansas  storm. 

Figure  9  shows  a  convective  cloud  turret  illuminated  entirely  by  lightning  taking 
place  within  the  cloud.  Figure  10  shows  several  lightning  channels  visible  in  the 
creases  of  the  convective  cloud.  Figure  11  shows  faint  lightning  channels  that 
apparently  terminate  in  clear  air.  Figure  13  shows  the  lightning  activity  taking 
place  within  the  cloud  when  the  measurements  shown  in  Figure  12  were  made. 


Optical  and  Electric  Field-Change  Data 

Data  acquired  during  the  Atlanta,  Georgia,  storm  on  5/11/81  are  of  uniformly 
better  quality  than  the  data  from  the  Arkansas  storm.  An  excellent  example  of  the 
optical,  E-field,  and  line-scan  camera  data  for  a  flash  to  ground  is  shown  in  Figure 
12(0),  and  in  the  expanded  plots  for  each  of  the  five  return-strokes  in  Figures  12(1) 
to  12(5).  It  should  be  noted  that  the  line-scan  spectrometer  is  a  charge-coupled 
device  that  accumulates  charge  until  it  is  read  out ,  the  video  signal  consisting  of  a 
contiguous  set  of  "frames"  that  recur  at  the  rate  of  1  every  5  msec.  The  signal  in 
a  particular  frame  is  due  to  the  charge  that  was  accumulated  by  the  device  during 
the  previous  5  msec.  Thus,  as  can  be  seen  in  the  figures,  the  spectra  usually 
follow  the  optical  lightning  event  by  as  much  as  5  msec.  It  should  also  be  noted 
that  the  amplitudes  of  the  optical  signals  of  the  individual  strokes  are  not  necessarily 
represented  faithfully  in  Figure  12(0).  The  problem  is  due  to  the  use  of  a  low 
sampling  rate  during  the  analog- to-digital  conversion  of  the  data  for  that  figure. 

The  true  relative  amplitudes  of  the  stroke  signals  can  be  inferred  from  the  expanded 
Figures  12(1)  to  12(5). 
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Figure  9.  Convective  cloud  turret  of  Georgia  storm  (approximately  11  km 
msl  altitude)  illuminated  by  lightning  within  the  cloud. 
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Figure  10.  Lightning  channels  visible  in  creases  of  cloud. 
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Figure  11.  Rather  faint  lightning  channels  apparently  terminating  in  clear  air. 
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Figure  12.  Optical,  electric  field,  and  spectral  data  from  a  flash  to  ground  over 
Atlanta,  Georgia  (0).  Signals  from  the  individual  strokes  are 
shown  expanded  in  (1)  through  (5). 
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dueed  the  data  shown  in  Figure  12* 


The  E-field  changes  produced  by  the  5  return-strokes  shown  in  Figure  12(0) 
are  each  accompanied  by  an  optical  signal,  and  as  in  Figure  5(0)  of  the  Arkansas 
storm,  there  is  little  or  no  detectable  optical  (or  electrical)  activity  that  occurs  in 
the  intervals  between  the  strokes.  The  first  stroke  of  the  flash  shows  evidence  of 
a  stepped-leader  in  both  the  E-field  and  the  optical  signals.  This  is  evident  in 
Figure  12(0),  where  the  leader  activity  can  be  seen  to  last  for  about  25  msec,  con¬ 
sistent  with  measured  stepped-leader  durations.  In  the  expanded  Figure  12(1),  the 
optical  signal  for  the  first  stroke  exhibits  only  a  simple  growth  and  decay  (even  on  a 
more  expanded  writeout).  The  first  stroke  is  not  the  strongest  of  the  five  strokes 
in  either  optical  or  electric  field-change  amplitude.  No  spectral  signals  are  above 
noise  in  Figure  12(1). 

In  Figure  12(2)  showing  the  second  stroke  we  see  an  excellent  example  of  the 
optical  representation  of  the  dart  leader  event.  There  is  excellent  correlation  with 
the  dart-leader  field  change,  and  the  light  signal  was  strong  enough  for  the  spec¬ 
trometer  to  register  a  spectrum  of  the  stroke  plus  leader  event. 

The  third  stroke  in  Figure  12(3)  produced  the  largest  field  change  and  the 
largest  light  output  of  all  five  strokes.  The  spectral  signal  is  very  strong,  showing 
essentially  only  the  two  dominant  lines  in  the  lightning  spectrum :  the  emission  lines 
of  singly  ionized  nitrogen  at  463.0  and  500.5  nm.  This  stroke  also  provides  an 
excellent  correlation  between  the  optical  dart  leader  and  the  dart  leader  field  change. 

The  fourth  stroke  [Fig.  12(4)]  is  very  much  like  the  second  stroke.  Again 
the  dart  leader  is  easily  distinguished  from  the  return-stroke  in  the  optical  signal. 

The  fifth  stroke,  shown  in  Figure  12(5),  is  small  in  amplitude  and,  coupled 
with  the  long  slow  field  change,  is  suggestive  of  a  long  continuing  current.  The 
oscillations  on  the  latter  part  of  the  field  change  are  spurious.  Despite  the  rather 
small  amplitude  of  the  discrete  light  pulses,  there  is  enough  integrated  light  intensity 
to  produce  a  signal  on  the  line-scan  spectrometer  trace,  a  fact  which  is  further 
evidence  for  the  existence  of  a  long  continuing  current. 

Additional  flashes,  both  intracloud  and  cloud- to- ground,  were  recorded  on  this 
flight.  They  are  essentially  similar  to  the  flashes  already  discussed  and  will  there¬ 
fore  not  be  presented. 


Spectral  Data 

Spectral  recordings  with  the  Fairchild  line  scan  camera  consistently  show  the 
two  strong  emission  features  recorded  in  the  photographic  data  and  reproduced  in 
Figure  8.  As  mentioned  in  the  previous  section,  these  features,  centered  at  463.0 
and  500.5  nm,  are  produced  by  emissions  from  singly  ionized  nitrogen.  No  lines 
were  observed  that  were  not  recorded  in  the  photographic  spectra  and  listed  in 
Table  1. 


DISCUSSION  AND  CONCLUSIONS 


With  few  exceptions,  the  photographs  that  indicate  the  presence  of  lightning 
within  the  cloud  show  hard  cauliflower-like  structures  characteristic  of  vigorous 
convection.  This  appears  to  offer  further  confirmation  of  Workman  and  Reynolds' 
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observation  [24]  that  strong  convective  activity  is  a  requirement  for  the  electrifica¬ 
tion  process  responsible  for  causing  lightning. 


In  some  cases,  such  as  those  shown  in  Figures  3  and  11,  it  is  clearly  evident N 
that  branches  of  the  lightning  channel  are  terminating  in  clear  air  above  some  por¬ 
tions  of  the  convective  cloud.  These  channels  are  quite  faint  in  comparison  with 
other  much  “brighter  onesT'fKat  bridge  creases  in  the  cloud  or  with  cloud-to-ground 
strokes  photographed  from  the  ground.  They  probably  transport  only  small  amounts 
of  charge.  Conceivably  there  are  other  very  weak  discharges  above  the  cloud  that 
are  too  weak  to  be  seen. 


It  is  impossible  to  tell  from  these  single  photographic  images  whether  the  dis¬ 
charge  extends  vertically  upward  above  the  cloud  or  whether  it  is  proceeding  either 
horizontally  or  even  downward  from  the  point  that  it  issued  from  the  cloud.  If 
second  photographs  of  these  channels  viewed  from  another  location  can  be  taken  froi 
another  airplane  or  from  the  gorund,  it  will  be  possible  to  determine  the  three- 
dimensional  structure  of  these  discharges. 

The  detailed  fine  structure  of  the  cloud  tops  that  can  be  seen  in  Figures  1  and 
9  is  not  unexpected.  It  has  long  been  evident  from  ground-based  observations  and 
in  particular  from  photographs  taken  from  high  altitude  aircraft  during  daylight  [25] 
that  the  penetrative  turrets  of  thunderclouds  have  a  complicated  cauliflower-like 
structure.  The  surprising  feature  .that  is  revealed  when  the  cloud  is  illuminated  by 
lightning  from  within  is  that  the  creases  and  folds  are  not  shallow  structures ,  but 
extend  downward  into  the  cloud,  possibly  as  much  .as. -several  kilometers."  Ubviously 
this  structure  indicates  the  presence  of  complicated,  organized  convective  air  motions 
They  probably  play  an  important  role  in  cloud  microphysical  processes  and  in 
thunderstorm  electrification . 


Studies  of  the  optical  pulses  and  the  electric  field  changes  accompanying  lightn¬ 
ing  as  viewed  from  above  the  storm  indicate  that  a  great  deal  of  specific  information 
regarding  lightning  events  (leaders,  return-strokes,  continuing  currents,  and 
spectra)  is  available  from  the  scattered  light  that  finds  its  way  out  of  the  top  of  the 
thundercloud. 


It  is  conceivable  in  this  small  sample  of  flashes  to  ground  that  the  optical 
sensor  was  viewing  the  "naked  channel"  below  the  cloud  as  well  as  the  scattered 
light  from  the  cloud  surface,  and  that  the  scattered  light  alone  would  not  have 
revealed  the  excellent  detail  in  the  dart  leader,  return-stroke  pairs.  To  show  that 
the  signals  presented  in  Figures  12(0)  through  12(5)  are  the  result  of  scattered  light 
alone  emerging  from  the  cloud  surfaces,  Figure  13,  which  was  taken  during  the  time 
of  the  flash  to  ground  identified  in  Figure  12,  is  given. 

The  situation  is  complicated  by  the  occurrence  of  another  discharge  at  the  time 
the  film  was  being  transported  after  the  event  of  Figure  12(0)  which  occurred  a 
03:32:32.  The  LED's  print  the  time  on  the  film  for  2  sec  prior  to  advancing  the 
film.  The  film  is  automatically  advanced  every  11  sec.  In  this  case,  the  frame  con¬ 
taining  the  photograph  of  interest  was  labelled  03:32:35,  indicating  that  the  film  was 
transported  into  position  at  03:32:26,  9  sec  before  the  LED's  were  activated.  This 
9  sec  period  contains  the  cloud-to-ground  flash  event.  An  examination  of  the  frame 
prior  to  the  one  of  interest  confirms  the  conclusion  that  the  displaced  channels 
present  in  the  photo  were  photographed  after  the  time  03:32:37,  when  the  film  was 
transported  into  the  next  frame.  Because  all  the  channels  visible  are  repetitive  and 
substantially  displaced  from  each  other,  they  must  have  all  occurred  during  the  150 
msec  interval  it  takes  to  transport  the  film .  The  displacement  is  not  due  to  the 
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motion  of  the  aircraft  which  flies  at  a  speed  of  200  m  sec  .  The  frame  area  on  the 
film  represents  a  field  of  view  of  about  10  x  10  km  at  cloud- top  height.  For  this 
discharge  (03:32:32),  which  lasted  for  less  than  300  msec,  the  maximum  displacement 
between  any  two  channel  images  would  have  to  have  been  less  than  60  m,  which 
amounts  to  a  barely  detectible  broadening  of  the  image. 

With  the  above  analysis,  the  existence  of  any  visible  lightning  channels  asso¬ 
ciated  with  the  discharge  at  03:32:32  can  be  ruled  out.  The  optical  signals  in 
Figures  12(0)  through  12(5)  were  therefore  the  result  of  scattered  light  emerging 
from  the  cloud  surfaces  and  not  associated  with  channels  in  clear  view  of  the  sensor. 

Although  there  is  time  dilation  produced  by  the  scattering  and  channel-lengthen 
ing  process,  it  is  clear  from  the  analysis  that  the  delays  and  superpositions  of  light 
signals  are  not  of  sufficient  duration  to  mask  the  identity  of  leaders  and  return 
strokes.  The  leader  luminosity  always  appears  to  be  less  than  that  of  the  return- 
stroke;  the  dart  leader  occurs  as  a  discrete  optical  event  that  precedes  the  return- 
stroke  by  1.5  msec  or  (usually)  less;  return-stroke  optical  events  in  a  flash  are 
usually  separated  by  periods  of  tens  of  milliseconds  with  practically  no  luminous 
events  between  them.  All  of  these  characteristics,  based  upon  optical  measurements 
alone,  can  be  used  to  write  a  pattern  recognition  algorithm  that  would  provide  a 
means  of  selecting  and  identifying  the  ground  flashes  in  a  storm,  and  to  sort  them 
out  from  the  intracloud  flashes.  [A  comparison  of  the  optical  signal  from  an  intra¬ 
cloud  flash  (Fig.  6)  with  that  of  a  cloud-to-ground  flash  [Figs.  5(0)  and  12(0)]  show 
a  striking  difference  between  the  flash  types  in  the  occurrence  frequency  of  pulses 
of  roughly  similar  amplitudes .  ]  A  reasonable  estimate  of  the  number  of  individual 
return- strokes  would  also  be  obtained. 

The  present  analysis  appears  also  to  have  identified  a  number  of  continuing 
current  strokes  through  the  use  of  a  light  integration  technique  (line-seal  spectrome¬ 
ter)  sampling  at  5  msec  intervals.  More  data  are  needed  here,  but  an  intensity  ratio 
between  the  two  lines  of  singly  ionized  nitrogen  is  clearly  evident.  Whether  this 
ratio  changes  between  intracloud  strokes,  high-current  return-strokes,  and  con¬ 
tinuing  current  strokes,  needs  to  be  studied.  If  it  does,  then  clearly  another 
method  exists  for  distinguishing  return-strokes  from  intracloud  flashes. 


SUMMARY 


Observations  of  nocturnal  thunderstorms  made  from  above  by  a  NASA  U-2 
airplane  disclose  features  of  cloud  structure  and  lightning  that  are  not  gene 
visible  from  the  ground.  Most,  if  not  all,  lightning  activity  seems  to  be  associated 
with  clouds  having  strongly  convective  cauliflower  tops.  When  the  lightning  occurs 
deep  within  the  cloud,  it  can  be  seen  that  the  folds  and  creases  of  the  cauliflower 
top  are  not  superficial  surface  characteristics,  but  extend  deep  into  the  cloud. 

In  both  of  the  storms  that  were  studied,  lightning  channels  are  visible  in  the 
clear  air  above  the  cloud  in  5  lo_J.O  percent  of  the  lightning  events.  T{iis  finding 
showS~1:hat  substances  known  tQ  be  produced  by  thunderstorm  electrical  discharges 
sucir  hs  "NO  and  0„  can  be  introduced  directly  into  the  stratosphere.  More  detailed 

measurements  will  be  required  to  determine  whether  or  not  the  quantities  introduced 
are  of  importance. 
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The  cause  and  nature  of  the  discharges  above  the  cloud  are  not  clear.  Possibly 
they  may  be  produced  by  accumulations  of  space  charge  in  the  clear  air  above  the 
cloud.  Alternatively  the  discharges  may  arise  solely  because  of  the  intense  electric 
fields  produced  by  charges  within  the  cloud.  In  the  latter  case  the  ions  introduced 
by  these  discharges  will  have  the  effect  of  increasing  the  electrical  conductivity  of 
the  air  above  the  cloud  and  will  increase  the  conduction  current  that  flows  from  the 
cloud  to  the  electrosphere. 

Lightning  spectra  observed  from  above  thunderstorms  do  not  appear  to  differ 
significantly  from  similar  observations  of  lightning  that  have  been  made  from  the 
ground  beneath  the  cloud.  More  quantitative  data  at  higher  resolution  may  eventually 
show  significant  spectral  differences  between  cloud-to-ground  and  intracloud  strokes. 

Our  findings  show  that  electric  field-change  data  taken  with  an  electric  field- 
change  meter  mounted  in  an  airplane  provide  data  on  lightning  discharges  from  above 
that  are  quite  similar  to  those  obtained  from  the  ground  in  the  past.  The  transients 
produced  by  the  photocell  optical  system  correlate  well  with  those  provided  by  the 
E  field-change  meter.  It  appears  that  even  though  E  field-change  measurements 
cannot  be  made  from  a  satellite,  the  optical  signals  from  dart  leaders,  from  return- 
strokes .  and  from  continuing  currents  are  recognizable,  and 'in  some  instances,  can 
be  used"t<rprovide  detailed  information  on  the  fine  structure  of  lightning.  In  par¬ 
ticular  ,  they  can  be  used  to  distinguish  between  cloud-to- ground  and  intracloud 
flashes. 

This  report  has  emphasized  the  "optical  alone"  aspect  of  sensing  lightning 
discharges  from  above  clouds  because  of  the  interest  in  satellite  measurements, 
especially  from  geosynchronous  altitude.  Clearly,  optical  measurements  from  a  single 
platform  are  uniquely  suited  to  provide  the  space  and  time  resolution  necessary  to 
view  an  Earth  hemisphere  on  which  hundreds  of  thunderstorms  may  be  active  simul¬ 
taneously.  Additional  scientific  objectives  have  been  discussed  in  Davis,  et  al.  [16], 
and  an  update  of  the  presently  used  instrumentation  package,  with  some  preliminary 
statistics  on  optical  pulse  frequency,  is  given  in  Christian,  et  al.  [26]. 
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